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INTRODUCTION 


The  following  report  represents  a  summary  of  the  work  that  has  been  done 
in  the  past  year  in  the  study  of  seismic  wave  scattering.  The  joint  CSDL/MIT 
efforts  in  this  study  have  been  separated  mto  three  areas  of  investigation.  Two 
areas  are  within  the  realm  of  ultrasonic  modeling.  First,  simple  models  are 
used  for  the  empirical  analysis  of  scattering  and  dispersion  Second,  structural 
features  found  in  the  earth  are  carefully  modeled,  making  it  possible  to  gen¬ 
erate  realistic  synthetic  seismograms.  The  third  area  of  interest  is  numerical 
modeling  where  important  contributions  have  been  made  through  finite 
difference  calculations. 

In  the  past  year,  the  work  on  these  three  fronts  has  begun  to  coalesce 
Insights  have  been  gained  concerning  the  propagation,  dispersion,  and  scatter¬ 
ing  characteristics  of  surface  waves  passing  through  a  medium  with  lateral  and 
surface  heterogeneities  The  research  has  also  received  a  clearly  defined  direc¬ 
tion.  In  the  future,  numerical  calculations  and  simple  lab  models  will  be  used  to 
quantitatively  describe  the  scattering  parameters  for  a  broad  class  of  different 
structural  and  topographic  geometries  In  conjuration  with  this,  detailed 
models  and  complicated  numerical  simulations  will  allow  geographic  localities 
to  be  modeled  in  the  lab.  This  will  have  important  consequences  concerning  the 
determination  of  regional  seismic  safety 

A  summary  of  past  research  in  wave  scattering,  using  analytic,  numerical, 
and  laboratory  methods,  is  given  below.  Various  analytic  treatments  of  lateral 
heterogeneities  have  been  given  by  Hudson  and  Knopoff  (1964),  Mai  and  Knopoff 
(1966),  McGarr  and  Alsop  (1967),  and  Gilbert  and  Knopoff  (1960) 

The  finite  difference  method  has  been  applied  to  scattering  by  Boore 
(1970),  Alterman  and  Loewenthal  (1972),  Munasinghe  and  Farnell  (1973),  Fuyuki 
and  Matsumoto  (1980),  and  Boore  et  a l  (1900) 
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Experimental  modeling  using  specific  geometries  can  be  found  in  Knopoff 
and  Gangi  (I960),  Martel  et  a L.  (1977),  and  Nathman  (1980). 

In  this  report,  the  contributions  from  experimental  and  numerical  work 
will  be  summarized  We  will  outline  the  work  that  has  been  accomplished  in 
each  of  the  three  areas  of  investigation  that  were  discussed  above  We  will  take 
a  particular  geometry  (the  case  of  a  rectangular  mountain'  or  valley  )  and 
show  how  it  is  treated  in  each  area  of  research.  Thus,  we  will  demonstrate  that 
the  results  from  each  area  complement  one  another. 


SIMPLE  LABORATORY  MODELS 

Piezoelectric  transducers  are  used  as  sources  and  receivers.  A  magnetos- 
trictive  material  such  as  nickel  is  used  as  the  medium.  The  case  of  steps  (Fig¬ 
ure  1)  has  received  careful  analysis.  The  seismogram  of  the  receiving  trans¬ 
ducer  can  be  digitized  for  subsequent  decomposition  of  the  waves.  For  exam¬ 
ple,  reflection  and  transmission  coefficients  can  be  calculated  making  it  possi¬ 
ble  to  compare  steps  of  different  height.  From  the  amplitude  and  phase  spectra 
we  can  identify  wave  composition  and  dispersion  characteristics.  Some  of  the 
results  can  be  found  in  Toksoz  et  a l.  (1981).  Similar  work  has  also  been  done 
by  Gangi  and  Wesson  (1978). 

Scattering  can  be  treated  quantitatively  with  this  model.  Receivers  can  be 
placed  as  shown  In  Figure  2  with  the  resulting  signals  shown  in  Figure  3.  The 
scattering  characteristics  of  different  shaped  obstacles  can  thus  be  catalogued. 

Varying  the  dimensions  of  rectangular  ’mountains'  is  possible  in  this  lab 
set-up.  The  results  are  shown  in  Figures  4  and  5.  It  is  observed  that  the 
transmitted  and  reflected  waves  are  sensitive  to  the  geometry  of  the  ’moun¬ 
tain’.  However,  detailed  study  of  the  differing  wave  components  is  not  possible 
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using  this  model  alone.  Comparing  the  reflections  from  rectangular  and  tri¬ 
angular  ’mountains’  (Figure  6),  we  can  deduce  that  a  larger  fraction  of  the 
energy  is  converted  into  body  waves  at  the  corner  for  the  rectangular  ’moun¬ 
tain*. 

ULTRASONIC  MODELING 

Much  of  the  work  at  CSDL  has  been  devoted  to  the  fabrication  of  functional 
transducers  and  multilayered  models.  A  detailed  summary  of  this  work  has 
been  given  in  the  CSDL  report  by  Chamuel  (1982)  which  is  included  in  this 
report  as  Appendix  1.  The  response  of  the  piezoelectric  transducers  that  have 
been  developed  is  shown  in  Figure  7.  An  important  characteristic  of  the  shear 
wave  transducer  is  shown  in  Figure  8:  the  capability  of  separating  transverse 
and  radial  components. 

The  seismic  properties  of  numerous  composite  materials  have  been  tested 
in  an  effort  to  suitably  model  sedimentary  layers.  The  interfaces  between  these 
materials  have  also  been  considered. 

The  transducers  and  modeling  materials  were  applied  to  the  three- 
dimensional  model  of  Dry  Lake  Valley,  Nevada  (shown  in  Figure  9).  This 
represents  an  extension  of  the  simple  rectangular  ’valley’.  The  model  was 
tested  with  and  without  sedimentary  layers,  thus  making  it  possible  to  distin¬ 
guish  between  the  dispersive  and  scattering  characteristics  of  layered  media 
and  topographical  features.  It  was  found  that  for  the  sedimentary-fUled  ’valley' 
the  transmission  of  Rayleigh  waves  was  greater  than  for  the  unfilled  ’valley’.  It 
is  also  useful  to  compare  Rayleigh  wave  transmission  through  a  flat  region  to 
that  across  an  empty  ’valley’.  The  ’valley’  has  the  effect  of  reducing  the  signal 
by  a  factor  of  10,  attenuating  higher  frequency  components,  and  creating  an 


earlier  body  wave  arrival. 

It  is  possible  to  model  the  rectangular  'mountain'  by  placing  two  solid 
layers  over  a  half-space.  In  figure  10  the  seismic  properties  of  the  model  are 
given  (with  the  solid  layer  parameters  given  in  the  first  two  rows).  The  seismo¬ 
grams  generated  by  displacing  the  receiver  away  from  the  layered  materials  are 
also  shown  in  Figure  10.  The  dispersive  characteristics  resulting  from  the  lay¬ 
ered  media  are  observed  in  these  traces.  These  traces  can  be  compared  to  the 
seismograms  in  Figure  4  for  the  case  of  homogeneous  layers.  Reflection  from 
layers  and  the  formation  of  new  modes  result  in  the  dominant  features  of  Fig¬ 
ure  10. 

NUMERICAL  MODELLING 

The  finite  difference  scheme  with  appropriate  boundary  conditions  has 
been  investigated.  Some  of  the  stability  problems  that  are  encountered  are 
described  in  Ilan  and  Lowenthal  (1976).  The  cases  of  steps  and  rectangular 
'mountains’  have  been  treated 

Figure  1 1  shows  the  input  surface  displacements.  The  Rayleigh  wave 
scattering  for  the  rectangular  mountain'  case  is  described.  The  geometry  is 
shown  in  Figure  12.  The  top  edge  is  a  free  surface  with  the  other  edge  modeled 
as  an  absorbing  boundary.  For  the  specifics  of  the  absorbing  boundary  condi¬ 
tions  see  Clayton  and  Engguist  (1977).  The  finite  difference  scheme  makes  it 
possible  to  monitor  the  wave  propagation  over  time.  In  Figures  13(a-i)  vectors 
representing  the  particle  displacement  are  plotted  at  each  grid  space  for  a  par¬ 
ticular  time  frame.  The  initial  waveform  is  shown  in  Figure  13a.  Many  insights 
into  wave  conversion  are  gained  from  this  simulation.  We  observe  the  reso¬ 
nance  of  the  ’mountain’.  Subsequently,  the  ’mountain’  acts  as  a  source  for 
body  waves.  The  seismograms  for  stations  on  the  incident  side  of  the  moun¬ 
tain’  are  shown  in  Figure  14  with  the  seismograms  from  the  transmitted  side 
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given  in  Figure  15.  Analysis  of  the  seismograms  and  ’snapshots'  of  Figure  13 
gives  us  information  on  the  wave  composition  in  the  near-  and  far-fleld 


The  following  observations  can  be  made  In  Figure  '3c,  the  Rayleigh  wave 
appears  to  track  along  the  free  surface.  The  first  cycle  of  the  Rayleigh  wave  has 
rotated.  In  Figure  13d  the  entire  ’mountain’  starts  to  resonate.  A  body  wave  is 
forming  at  the  right  corner  of  the  mountain’.  The  surface  wave  itself  appears 
to  float  across  the  inside  of  the  ’mountain’  just  making  contact  with  the  other 
side  It  is  unclear  if  this  energy  would  be  considered  a  body  wave  or  a  surface 
wave  when  it  crosses  the  interior  of  the  ’mountain’  The  ’mountain'  continues  to 
resonate  as  the  energy  trapped  in  the  ’mountain’  reverberates.  In  Figure  13f,  a 
body  wave  has  formed  at  the  left  corner  of  the  ’mountain'  while  the  transmitted 
surface  continues  to  move  away  from  the  mountain'.  The  ’mountain’  continues 
to  resonate  In  Figures  13(g-i),  the  Rayleigh  wave  leaves  the  grid  and  the 
'mountain'  finally  damps  out  by  radiating  energy  as  body  waves  to  the  lower 
media 

Figure  14  shows  seismograms  1-3  from  the  incident  side  of  the  'mountain 
as  found  in  Figure  1  3.  Seismogram  3  is  closest  to  the  mountain’  Complicated 
reflected  phases  can  be  seen  along  with  the  incident  wave  Figure  *4  also  shows 
seismograms  4-8  on  the  transmitted  side  The  amplitudes  gradually  increase  to 
a  stationary  value  away  from  the  'mountain'  Also,  the  pulse  compresses  away 
from  the  step.  In  seismograms  7-8,  a  body  wave  is  seen  to  form  m  front  of  the 
transmitted  surface  wave. 


C  'NCU  Wn 

Ultrasonic  memods  applied  to  scale  models  of  realistic  earth  structures 
and  simple  representative  structures  offer  a  powerful  means  of  investigating 
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scattering  phenomena.  In  addition,  the  finite  difference  method  provides  a  use- 
ful  new  approach  to  the  interpretation  of  seismic  wave  scattering.  For  example, 
the  interaction  of  a  Rayleigh  wave  with  a  simple  model  'mountain'  is  demon¬ 
strated  with  this  method  to  be  a  complex  process  in  which  the  Rayleigh  wave 
excites  oscillations  of  the  'mountain*  as  a  unit.  These  oscillations  decay  by  re¬ 
radiation  of  body  wave  energy  into  the  medium. 

We  now  have  the  means  to  interpret  ultrasonic  seismograms  using  tech¬ 
niques  similar  to  those  used  in  terrestrial-scale  seismology.  The  expected 
result  is  that  laboratory  modeling  will  receive  a  secure  quantitative  foundation, 
enabling  correspondence  between  laboratory  models  and  earth  structure  to  be 
made  with  greater  confidence  in  the  future. 
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Figure  1.  Laboratory  set-up  for  two-dimensional  step 
model  . 


Figure  2.  Receiver  positions  used  to  measure  azimuthal 
scattering  of  step. 
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:  Waveforms  representing  the  input  pulse  for  the  finite 
difference  model. 


Fig.  2 -  Area  of  interest  for  finite  difference  model. 
Receiver  locations  are  also  shown. 
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Fig.  13:  Time  sequence  depicting  the  scattering  of  a  surface 
wave  for  the  rectangular  ‘mountain1  case. 
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THREE  DIMENSIONAL  SEISMIC  ULTRASONIC  MODELING 
PART  I  :  TECHNIQUES  AND  EXPERIMENTAL  RESULTS 

Jacques  R.  Chamuel 

The  Charles  Stark  Draper  Laboratory,  Inc. 


INTRODUCTION 

The  seismic  vulnerability  of  ground  facilities  such  as  nuclear  power  plants 
and  military  installations  depends  on  ground  motion  duration  and  amplitudes.  It  is 
essential  to  understand  and  be  able  to  predict  these  ground  motions  generated  by 
earthquakes  and  explosions. 

Real  three  dimensional  seismic  problems  are  practically  beyond  the  state 
of  the  art  of  theoretical  formulation  and  numerical  analysis.  In  qeneral,  seismic 
problems  are  too  costly  to  be  solved  with  reasonable  speed  and  accuracy  with  nu¬ 
merical  techniques.  Furthermore,  empirical  determination  by  field  measurements 
of  the  effects  of  the  earth  structure  on  the  ground  motion  caused  by  large  events 
(nuclear  explosions)  are  not  practical. 

In  1979,  Jacques  Chamuel  (CSDL)  demonstrated  the  potential  feasibility 
of  three  dimensional  seismic  ultrasonic  modeling  by  using  a  carved  limestone  block 
filled  with  wax  composite  materials.  The  limestone  block  was  a  free  sample 
provided  by  C.S.  Precourt  &  Sons,  Co.,  Sudbury.  Prof.  Nafi  Toksoz  (MIT)  ) 

and  Prof.  Howard  Schif f le  .t v 3 *  (Long  Beach  State  College,  CA)  participated  as 
consultants  to  implement  Chamuel1 s  proposed  seismic  ultrasonic  modeling  techniques 
to  develop  a  quantitative  model  of  Dry  Lake  Valley,  Nevada. 


THREE-DIMENSIONAL  SEISMIC  ULTRASONIC  MODELING  TECHNIQUES 
A.  Modeling  Materials 

In  past  work,  the  Charles  Stark  Draper  Laboratory,  Inc.  developed  a  carved  scaled 
model  of  Dry  Lake  Valley  region,  Nevada.  A  photograph  of  the  unfilled  model  (scale 
1:250,000)  is  shown  in  figure  1.  The  actual  dimensions  of  the  model  are  38x28x21  cm. 
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Bedford  limestone  was  chosen  to  model  the  rocks.  Under  the  current  AFGL  contract, 
CSDL  has  developed  necessary  modeling  materials  to  fill  the  valleys  of  the 
carved  model.  The  preliminary  materials  developed  under  the  BMC  task  were  not 
suitable  for  modeling  the  top  earth  layers  having  a  Rayleigh  velocity  between 
450  and  700  m/sec. 

Our  approach  to  the  modeling  materials  problem  has  been  tc  look  for 
materials  that  can  provide  us  with  a  wide  ranges  of  velocities  ind  .iensities, 
and  at  the  same  time  are  easy  to  handle  and  replace.  Epoxies  in  general 
create  toxic  fumes  and  cannot  be  removed  easily  from  the  carved  model  for 
modification  or  replacement.  We  decided  to  pursue  the  wax-base  modeling 
materials  since  they  are  in  many  respects  practical  for  preliminary  three-dimen¬ 
sional  modeling. 

Figure  2  shows  the  simplified  layered  structure  of  Dry  Lake  ’/alley  and  tne 
preliminary  modeling  materials  used.  Table  1  summarizes  the  measured  properties 
of  the  seismic  ultrasonic  modeling  materials.  In  our  simplified  model,  the 
valley  was  filled  with  two  parallel  layers.  The  top  layer  was  modeled  using 
paraffin  wax  mixed  with  vaseline  and  aluminum  powder.  The  middle  layer 
material  consisted  of  Camauba  wax  and  silica  powder.  A  2%  (by  weight)  amount 
of  beeswax  was  added  to  the  middle  layer  material  to  reduce  the  shear  wave 
velocity.  As  can  be  seen  from  Figure  2,  we  were  able  to  match  the  modeling 
material  properties  to  the  actual  earth  properties.  The  top  layer  thickness 
was  1.7  mm,  and  the  middle  layer  filled  the  gap  between  the  bottom  surface 
of  the  top  layer  and  the  limestone  block.  Away  from  the  deepest  spot  of  the 
valley,  the  thickness  of  the  middle  layer  was  2.3  mm. 

Although  the  described  modeling  materials  have  many  advantages,  including 
processing  in  a  relatively  short  period  of  time,  the  elastic  properties  of 
these  composite  waxes  are  highly  dependent  on  temperature.  We  made  preliminary 
tests  to  determine  their  velocities  temperature  dependence.  The  Rayleigh  wave 
velocity  of  surface  waves  propagating  along  the  edge  of  modeling  materials 
plates  was  monitored  at  different  temperatures.  For  the  soft  top  layer  material 
{vaseline-paraffin-aluminum  compound)  Rayleigh  velocities  were  determined  to 
be  690  and  463  m/sec  at  65  and  85  degrees  Fahrenheit,  respectively.  The  middle 
layer  material  is  much  less  sensitive  to  room  temperature  variations.  The 
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Figure  2  .  Preliminary  modeling  materials  used  to  match  seismic  properties  of  Dry  Lake  Valley. 
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The  effects  of  temperature  on  a  layered  simple  model  are  illustrated  in 
Figure  3.  The  model  is  made  of  two  parallel  layers  on  a  large  limestone  block . 
The  top  layer  was  69  mil  thick,  and  the  middle  layer  was  91  mil  thick.  Again 
we  used  the  same  materials  described  above  for  the  layers.  The  effects  of 
raising  the  temperature  of  the  model  on  the  dispersion  and  attenuation  ot  the 
received  signal  may  be  observed  from  the  wave  trains  of  Figure  3.  Raising  the 
temperature  softens  the  top  layer  material  and  decreases  it  velocities.  The 
time  of  arrival  of  the  high  frequency  components  would  vary  accordingly.  It 
is  clear  that  a  temperature  controlled  environment  is  needed  to  achieve  accurate 
modeling  with  the  soft  wax  modelling  materials. 


We  are  currently  evaluating  alternative  modeling  materials  for  the  layered 
valleys.  We  were  able  to  reduce  the  Rayleigh  wave  velocity  of  dtyeast  epoxy 
resins  by  adding  plasticizers. 
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Figure  3. 


Effects  of  model  temperature  on  Rayleigh  wave 
propagation  in  a  two- layer  over  half -space  model. 
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B.  Fabrication  of  Multilayered  Models 

In  order  to  fabricate  multilayer  models,  we  used  ultrasonic  waves  to  determine 
the  solidification  temperatures  of  the  wax-base  modeling  materials.  Ultrasonic 
wave  reflection  information  obtained  from  waves  travelling  in  a  magnetrostrictive 
rod,  immersed  in  the  modeling  material,  was  plotted  as  the  modeling  material 
temperature  decreased.  The  resulting  plots  of  the  middle  layer  and  the  top 
layer  materials  are  shown  in  Figure  4.  As  the  waxes  solidified,  the  wave  ref¬ 
lections  decreased,  since  the  solidified  wax  attenuated  the  travelling  waves. 

Both  temperature  and  reflection  measurements  were  recorded  as  the  material  was 
allowed  to  cool.  The  solidification  temperature  of  the  top  and  bottom  layer 
materials  were  95  and  175  degrees  Fahrenheit.  The  mixtures  remained  liquid 
down  to  162  and  220  degrees  respectively.  We  could  safely  fabricate  a  double 
layer  model  without  disturbing  the  bottom  layer  by  remelting. 


C.  Transducers 

Two  compressional  and  one  shear  wave  piezoelectric  transducers  were 
developed  during  the  course  of  the  study.  The  piezoelectric  elements  were 
P2T-5  having  a  4.7mm  diameter  and  a  2  M llz  equivalent  thickness.  Lead  was  used 
to  damp  the  reflections  of  the  elastic  waves  from  the  boundary  of  the  piezo¬ 
electric  elements  opposite  the  sensing  surface.  A  stainless  steel  cylinder  was 
machined  to  house  the  PZT  transducers .  A  photograph  of  the  three  transducers 
is  shown  in  Figure  1.  The  sensing  surface  was  covered  with  an  epoxy  composite 
material  to  protect  the  transducer.  The  diameter  of  the  sensing  area  was 
reduced  to  about  3  mm  by  machining  the  tip  of  the  epoxy.  The  response  of  the 
two  compressional  transducers,  when  placed  5  cm  apart  along  the  edge  of  a  1  mm 
thick  aluminum  sheet,  can  be  seen  in  Figure  2  (time  scale  5  usec/div) .  The 
compressional  wave  arrived  at  11  usee  and  the  Rayleigh  wave  at  about  18  usee. 

The  directivity  of  the  shear  transducer  was  tested  by  placing  the  shear 
transducer  perpendicular  to  the  plane  of  the  sheet  and  orienting  it  first 
radially,  then  transversely  to  the  wave  propagation  direction.  The  source  was 
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Figure  4.  ULTRASONIC  MONITORING  OF  SOLIDIFICATION 
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a  compressional  transducer  placed  at  the  edge  of  the  sheet.  The  relative 
transverse  and  radial  directional  response  of  the  shear  transducer  is  shown  in 
Figure  7  .  These  transducers  were  used  to  study  the  Rayleigh  wave  scattering 
from  two  and  three  dimensional  models  as  described  in  the  report. 


RAYLEIGH  WAVE  SCATTERING  FROM  VERTICAL  DISCONTINUITIES 

Rayleigh  waves  propagating  on  the  surface  of  the  earth  are  greatly  affected 
by  the  presence  of  topographical  discontinuities.  As  a  follow  on  to  our  studies 
on  the  effects  of  topography  on  Rayleigh  wave  scattering,  we  measured  from  the 
three  dimensional  Dry  Lake  Valley  model  the  transmitted  wave  propagating  across 
the  valley  without  sediment.  Figure  8  shows  the  outline  of  the  valley  drawn 
at  its  center  (scale  1:1).  The  received  signal,  Figure  8(b),  is  compared  with 
another  signal  received  from  a  half-space,  Figure  8(a)  keeping  the  same  horizontal 
distance  (4  cm)  between  the  source  and  receiver.  The  reflections  and  scattering 
from  the  valley  decreased  the  Rayleigh  wave  by  almost  a  factor  of  ten.  From 
Figure  8,  one  can  notice  the  reduction  of  the  high  frequency  components ,  the 
delay  in  the  Rayleigh  wave  arrival,  and  the  appearance  of  an  early  arrival  wave 
caused  by  conversion  of  body  waves  into  surface  waves. 

Figure  9  shows  the  transmitted  waves  obtained  from  a  two  dimensional  limestone 
model  consisting  of  a  simple  rectangular  cut  having  the  same  width  as  Figure  8 
Valley,  and  a  depth  of  7.5  mm  (simplified  model  of  Valley).  The  transmitted 
signal  also  consists  of  an  early  arrival  wave  followed  by  the  Rayleigh  Wave. 

A  comparison  of  Figure  9  with  Figure  8(b)  indicates  that  the  lack  of  the 
multiple  steps  in  the  approximate  model  resulted  in  a  larger  transmission  of 
the  high  frequency  converted  body  waves,  and  at  the  same  time  resulted  in  a 
reduction  of  the  high  frequency  components  of  the  Rayleigh  Wave.  It  would  be 
interesting  to  compare  these  experimental  results  with  a  theoretical  model  based 
on  the  superposition  of  steps. 
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Figure  9.  Rayleigh  wave  transmission  across  a  simplified  2-D  model 
of  Dry  Lake  Valley  center  (no  sediments) 


RAYLEIGH  WAVE  DISPERSION  FROM  DOUBLE  SOLID  LAYERS  OVER  HALF-SPACE 

A  double  layer  model  was  fabricated  and  tested  to  evaluate  our  seismic 
ultrasonic  modeling  technique.  A  2.3  mm  thick  layer  composed  of  the  camauba 
wax  mixture  was  deposited  on  a  limestone  block.  A  second  layer  was  added  on 
top  of  the  camauba  mixture  layer.  The  top  layer  was  1.6  mm  thick  and  consisted 
of  the  vaseline  paraffin  mixture.  Compressional  piezoelectric  transducers 
were  used  to  generate  and  detect  the  elastic  waves  in  the  double  layer  model. 
Initially,  the  receiver  was  1.75  cm  away  from  the  source,  then  it  was  displaced 
in  0.5  cm  intervals.  The  received  signal  was  photographically  recorded  from 
the  oscilloscope.  Figure  10  shows  the  family  of  signals  recorded  from  the 
experimental  model.  The  elastic  parameters  of  the  layers  are  given  in  Figure  10. 
The  group  velocity  experimental  points  were  calculated  and  plotted  in  Figure  11 
as  function  of  frequency.  The  theoretical  group  and  phase  velocities  will  be 
included  in  the  MIT  final  report.  Our  preliminary  analysis  indicates  that  we 
have  a  satisfactory  modeling  technique. 
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RAYLEIGH  WAVE  SCATTERING  from  FILLED  THREE  DIMENSIONAL  DRV  LAKE  VALLEY  MODEL 


The  unfilled  Dry  Lake  Valley  model  (Figure  1)  was  filled  with  two  parallel 
layers  of  sediments.  The  top  layer  is  69  mil  thick  composed  of  vaseline,  paraffin 
wax,  and  aluminum  powder.  The  second  layer  his  different  dimensions  depending 
on  the  depths  of  the  valley.  The  average  thickness  of  the  second  layer  is  about 
2.3  mm.  We  used  carnauba  wax,  beeswax,  and  silica  powder  to  model  the  second 
layer. 

The  elastic  properties  of  these  modeling  materials  are  identical  to  the 
double  layer  over  half  space  model  described  in  the  previous  section. 

A  photograph  of  the  filled  three  dimensional  model  with  the  transducers 
is  shown  in  Figure  12. 

Preliminary  Rayleigh  wave  measurements  have  been  made  using  this  model  to 
determine  the  effects  of  topography  on  Rayleigh  wave  scattering  and  dispersion. 

Figure  13  shows  the  relative  signals  obtained  as  the  Rayleigh  wave  path  crossed 
a  flat  region  (R^)  or  a  valley,  (R^  and  R^) .  The  time  scale  setting  was  10  usec/div . 
Transmission  across  the  valley  crested  dispersion,  made  conversion,  and  attenua¬ 
tions  . 

We  noticed  that  filling  the  valleys  increased  the  ground  motion  duration 
and  amplitude.  Figure  14  shows  the  effects  of  filling  a  valley  from  a  two 
dimensional  limestone  model.  Figure  15(a)  shows  the  effects  of  the  sediments 
on  the  transmission  across  the  three  dimensional  Dry  Lake  Valley.  The  received 
signal  is  much  larger  than  Figure  8(b)  where  the  valley  was  empty  (absolute 
amplitude  information  is  questionable) .  We  damped  the  Rayleigh  wave  of  Figure 
15(a)  by  pressing  by  hand  on  the  sediment  surface.  Figure  15(b)  shows  the 
received  signal  with  the  damped  Rayleigh  wave.  The  transmission  across  the 
filled  valley  in  Figure  13  (R^)  is  almost  1/4  of  Figure  13  (R^)  signal,  however, 
in  the  empty  valley  of  Figure  8(b)  the  ratio : transmission  across  a  valley/trans¬ 
mission  across  a  flat  region  was  about  1/10  in  our  model. 


10  usec/div. 
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Figure  14-  Rayleigh  wave  transmission  across  unfilled  (top)  and  filled 
(bottom)  two  dimensional  limestone  model  of  Dry  Lake  Valley 
center  line.  Used  paraffin  and  aluminum  as  sediment. 


Figure  15.  Rayleigh  wave  transmission  across  the  width  of 
Dry  Lake  Valley  model  with  transducers  outside 
of  valley  (  x=4cm) , 
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A  set  of  measurements  were  taken  from  the  filled  Dry  Lake  Valley  to  study 
Rayleigh  wave  dispersion.  The  source  was  located  in  the  valley  near  its  north 
end.  The  receiver  was  displaced  1  cm  at  a  time  towards  the  south.  The  initial 
distance  between  the  transducers  was  2  cm.  Figure  16  shows  the  dispersion 
results.  The  time  scale  was  10  usec/div.  At  8  cm,  an  additional  waveform 
with  a  20  usec/div  setting  was  recorded.  The  high  frequency  signals  observed 
at  about  110  usee  are  the  reflections  of  the  P  waves  from  the  bottom  surface 
of  the  large  limestone  block. 

We  used  the  shear  wave  transducer  to  measure  the  radial  and  transverse 
components  of  the  dispersed  waves.  Figure  (17)  was  taken  from  the  Dry  Lake 
Valley  model.  The  rejection  of  the  Rayleigh  Wave  is  observed  when  the  shear 
receiver  transducer  was  oriented  transversely. 

RECIPROCITY 

Rayleigh  waves  reflection  coefficient  is  much  larger  for  a  downstep 
topographical  surface  discontinuity  than  for  an  upstep  case.  Intuitively, 
one  may  think  that  interchanging  source  and  receiver  positions  would  give 
different  transmission  for  an  upstep  versus  a  downstep  case.  The  author  was 
curious  to  test  the  reciprocity  theorem  on  the  three  dimensional  model  in  view 
of  conversion  of  surface  to  body  waves  and  body  waves  to  surface  waves.  A 
source  was  placed  at  "SM  on  the  surface  of  the  filled  valley  as  shown  in 
Figure  19.  A  receiver  was  placed  at  "RM .  '»  recei  ved  signal  is  displayed 

in  Figure  20(a).  The  electrical  leads  of  the  transducers  were  interchanged 
so  the  original  source  is  now  used  as  a  receiver.  The  received  signal  for 
the  interchanged  case  is  shown  in  Figure  20(b).  The  same  experiment  was 
repeated  without  moving  the  oscilloscope  trace.  Figure  20(c)  shows  the 
superposition  of  the  two  signals.  A  close  look  at  Figure  20  signals  reveals 
that  Figure  20(b)  is  slightly  smaller  than  20(a).  This  is  due  to  the  different 
electrical  impedances  of  the  transducers.  The  experiment  was  repeated  on  a 
half-space  (limestone  block)  giving  the  same  ratio.  From  these  result,  we 
conclude  that  the  Reciprocity  theorem  holds  for  conversion  of  surface  and 
body  waves. 

In  the  real  earth,  the  presence  of  nonlinearities  would  tend  to  modi  tv 


the  problem. 
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Figure  17.  Radial  and  transverse  components  measured  from 

three-dimensional  model  with  source  and  receiver 
located  on  Dry  Lake  Valley  (x=5cm) • 
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Figure  19  .  Reciprocity  test  results:  (a)  source-receiver, 
(b)  receiver-source  ,  (c)  source-receiver 
superimposed  on  receiver-source  .  Data  taken 
from  3-D  model  of  Dry  lake  Valley.  Transducers 
locations  shown  in  figure  18. 


SUMMARY  AND  CONCLUSIONS 

The  Charles  Stark  Draper  Laboratory  has  continued  the  development  of  three- 
dimensional  seismic  ultrasonic  modeling  techniques  for  Air  Force  applications. 

New  broadband  piezoelectric  compressional  and  shear  transducers  were  developed. 

New  modeling  materials  suitable  for  modeling  top  earth  layers  were  introduced. 

A  preliminary  three-dimensional  model  of  Dry  Lake  Valley,  Nevada  region  with 
sediment  filled  valleys  was  constructed.  Experimental  results  showing  the  effects 
of  layers  and  topographical  features  on  Rayleigh  wave  attenuation,  dispersion, 
and  scattering  were  presented.  We  planned  and  conducted  experiments  demonstrating 
receiving  identical  signals  when  source  and  receiver  locations  are  interchanged 
fox  Rayleigh  wave  transmission  across  complex  topographical  features. 

The  main  conclusions  to  be  drawn  from  the  CSDL  efforts  are: 

1.  Filling  the  model  valleys  with  sediments  increased  the  Rayleigh  wave  transmission 
across  the  valley. 
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2.  The  wax-base  modeling  materials  are  suitable  for  fabrication  of  multilayered 
models.  Temperature  controlled  models  would  be  required  for  quantitative  studies. 

3.  Rayleigh  wave  transmission  across  an  empty  valley  (Dry  Lake  Valley)  decreased 

the  signal  by  almost  a  factor  of  ten,  reduced  the  high  frequency  content,  and 

created  an  early  arrival  wave  caused  by  body  waves. 

4.  The  presence  of  multiple  steps  in  the  valley  underneath  the  sediments  drastically 
reduces  the  high  frequency  content  of  the  early  arrival  transmitted  wave. 

5  deling  the  valley  with  one  major  slot  having  vertical  boundaries  (no  sediments) 
increased  the  transmission  of  the  high  frequency  converted  body  waves,  and  at 
the  same  time  resulted  in  a  reduction  of  the  high  frequency  components  of  the 
transmitted  Rayleigh  wave. 

6.  We  discovered  that  reciprocity  holds  for  conversion  of  body  waves  to  surface 

waves  and  also  for  conversion  of  surface  waves  to  body  waves  resulting  from 

Rayleigh  wave  transmission  across  topographical  features. 
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